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Abstract: We analyse several implications of lensing by the regular component of 



O 
O 
O 

p i' the galactic magnetic field upon the observed properties of ultra high energy cosmic 

Q ■ rays. Magnetic fields deflect cosmic ray trajectories, causing flux (de) magnificat ion, 

+i ■ formation of multiple images of a single source, and time delays. We derive the 

• energy dependence of these effects near the caustics at which the flux amplification 

of a point source diverges. We show that the large magnification of images around 
^ . caustics leads to an amplification bias, which can make them dominate the flux in 

■ some energy ranges. We argue that clustering in the arrival directions of UHECRs of 

comparable energy may be due to magnetic lensing around caustics. We show that 
magnetic lensing can also significantly alter the observed composition of cosmic rays 
at the highest energies. We also show that the time delay between events from a 
single image may monotonically decrease with decreasing energy in the neighborhood 
of a caustic, opposite to its behaviour in normal regions. 
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1. Introduction 

Galactic and intergalactic magnetic fields deflect extragalactic charged cosmic ray 
trajectories in their journey from their sources to the Earth, causing several effects 
upon the observed properties of ultra high energy cosmic rays (UHECRs). For a 
review see for instance [[l|. 

In a previous paper (in what follows paper 1), we have shown that the regular 
component of the galactic magnetic field acts as a giant lens upon charged CRs, and 
this can sizeably amplify (or demagnify) the flux arriving from any given source, 
modifying its spectrum. The magnification of the CR flux by the galactic magnetic 
field becomes divergent for directions along critical curves in the sky seen from the 
Earth, corresponding to caustic curves in the "source plane" , i.e. in the corresponding 
directions outside the Galaxy. The location of the caustics move with energy and, 
as a caustic crosses a given source direction, pairs of additional images of the source 
appear or disappear. Multiple image formation is a rather common phenomenon, at 
least within the galactic magnetic field models considered in 1. Indeed, the caustics 
sweep a rather significant fraction of the sky as the ratio E/Z between the CR energy 
and charge steps down to values of the order of a few EeV (1 EeV = 10^^ eV), before 
the drift and diffusive regimes turn on. At values of E/Z larger than around 50 EeV 
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caustics are present but sweep out a relatively small fraction of the sky. The effects 
under discussion are thus relevant even for the highest energy events so far detected 
if the CRs have a component which is not light. 

In the present paper we analyse in detail the energy dependence of the magnifi- 
cation of image pairs near a caustic, and discuss several implications of the existence 
of caustics upon the observed properties of UHECRs. 

In section ^ and in appendix A we show that if a source lies along a caustic 
at energy Eq, the magnification factor fi of each image in the pair that becomes 
visible at energies below Eq diverges as yU ~ A/ y^l — E /Eq. A is a dimensionless 
constant, fixed by the structure of the magnetic field along the CR trajectory to 
the source. It must be determined numerically, and we do so for some examples of 
source locations. We discuss several implications of this result, such as the expected 
enhancement in the observed event rate at energies near a caustic, which can lead 
to an amplification bias, making it more likely to detect images of UHECR sources 
that lie along caustics than sources in ordinary regions. In appendix B we provide a 
simple proof of the relation between the magnification of the three images appearing 
when there are two nearby folds. In section |^ we show that the angular separation 
between images in a pair increases near the caustic as /S.6 oc ^1 — E/ Eq. We analyse 
the expected event rate from the original source and its multiple images not only as 
a function of energy but also in terms of the observed arrival direction. We show 
that magnetic lensing near caustics is a source of clustering in the arrival directions 
of events with comparable energy. In section ^ we stress the fact that at fixed 
energy the magnetic lensing effects depend upon the CR electric charge, and show 
that this dependence can significantly alter the observed CR composition at the 
highest energies. In section || we discuss features of the time delay between events 
from multiple images, relevant in the case of bursting or highly variable sources of 
UHECRs P]. We show that the delay between the arrival of events of equal energy 
from the two images in a pair increases as At oc ^1 — E/Eq around a caustic. We 
also show that the time delay between events at different energy from one image in a 
pair may monotonically decrease with decreasing energy near the caustic, opposite to 
the behaviour of time delays in normal regions. Section |^ rounds up our conclusions. 

The galactic magnetic field model that will be used throughout this paper to 
illustrate magnetic lensing effects is the bisymmetric spiral model with even symme- 
try (BSS-S) described in paper I, which is a smoothed version of one of the models 
used in and |^ to study the effects of CR deflections in our galaxy. We refer the 
reader to paper I for details about the field configuration and about the numerical 
methods implemented to determine CR trajectories and flux magnifications. 

2. Flux enhancement by magnetic lensing near critical points 

As shown in paper I, the galactic magnetic field can act as a giant lens that amplifies 
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or demagnifies extragalactic sources of UHECRs, much alike the gravitational lensing 
effect upon distant quasars by intervening matter along the line of sight 0. CRs 
from a distant extragalactic source that enter the galactic halo from a direction b) h 
(in galactic coordinates) are deflected by the magnetic field and thus observed on 
Earth as if coming from a different direction {i,b)E, and their flux is amplified (or 
demagnified) by a factor /i. The magnitude of the effect depends upon the direction 
of observation and upon the ratio E/Z between energy and charge of the CR. The 
effect is most dramatic at the critical curves of the lens mapping, the directions along 
which the observed flux diverges for a fixed value of E/Z. The corresponding lines 
in the source coordinates (the direction from which CRs enter the galactic halo) are 
the caustics of the lens mapping. The location of the caustics changes with energy, 
for fixed Z. 

If a source position lies along a caustic of the magnetic lens mapping at energy 
Eq, a pair of images of the source either becomes visible or dissapears at energies 
below Eo. Their magnification diverges at E = Eq. This behaviour was illustrated in 
paper I for sources at galactic coordinates (£, b) = (282.5°, 74.4°) (M87 in the Virgo 
cluster) and {i,b) = (320°, —30°) (visible from the southern hemisphere). 

In appendix A we give a geometrical interpretation of the magnification of a 
pair of images near a caustic. Similar to the gravitational lens case ||, ||, the 
magnification diverges as 1/y/x, where x measures the distance of the source to the 
caustic. In the magnetic lensing case the location of the caustics "move" with energy, 
and thus the magnification of a pair of images diverges as the energy approaches the 
energy at which the source lies along the caustic as 1/ ^/Eq — E. In appendix A we 
also find the first two corrections to this leading order behaviour, and show that the 
magnification of the two images behaves as: 

f^iiE)^ , ^ ±B + aJl-E/Eo . (2.1) 
^l^E/E^ 

The dimensionless coefficient A is the same for the two images. The constant term 
has the same value with opposite sign for each image. The third term has a 
different coefficient Ci for each image due to the truncation of the expansion up to 
this order. The value of the coefficients A, B and Ci are fixed by the properties of 
the magnetic field along the CR trajectories. We determine them through a ffi to the 
numerical output for the energy dependence of the amplification, evaluated through 
the method described in paper I. 

Figure ^ displays the numerical result for the amplification near the caustics 
along with the ffi to the analytic expression in eq. ( p.l|) , for the same examples of 
source locations as in paper I. The top panel corresponds to M87 and the bottom 
panel to the source in the southern hemisphere. The figures display the magnification 
of the principal image (the one that is also visible at the highest energies) and of the 
images A and i?, visible at energies below the caustic only. The analytic expression 
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( |2.1|) fits with very high accuracy the numerical result for the magnification of the 
secondary images near the caustic, at least down to energies 10% below the energy 
Eo of the caustic. The fit to M87 determines Eq/Z = 20.41 EeV, A = 1.3, B = 4.0, 
Ca = —3.1, Cb = 5.0. The fit to the southern source fixes Eq/Z = 15.425 EeV, 
A = 0.44, B = 0.37, Ca = -0.19, Cb = -0.28. 
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Figure 1: Numerical result (points) and analytic fit to eq. (|2.l| ) (solid lines) for the 
amplification near a caustic of a pair of images of a point source. The source in the 
top panel is at {i,b) = (282.5°, 74.4°) (M87) and the source in the bottom panel is at 
{i,b) = (320°,— 30°). The numerical result for the amplification of each principal image 
(diamonds) is also shown. 

The divergence of the flux magnification of a CR source ai E = Eq is softened 
down to finite values when integrated across an extended source. This still allows 
for extremely large magnifications. The limiting factor to the maximum attainable 
magnification in a realistic situation is not the extended nature of the sources, but 
the fact that the divergence in the magnification, even in the case of a point source, 
arises only at a fixed energy Eq. Since realistic sources are not monoenergetic, the 
integrated fiux of a magnified source around Eq is thus always finite, even in the 
point source approximation. 
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Large magnification of image pairs around caustics leads to a significant enhance- 
ment of tlie detection probability of a source in a flux-limited sample. In gravitational 



lensing this effect is termed "amplification bias"|TO|, and may be responsible for the 
observation of quasars that would otherwise be too dim to detect. Let us now consider 
the strength of this effect due to galactic magnetic lensing of UHECRs. Consider a 
differential flux of CRs injected in the galactic halo at energies beyond the "ankle" 
given by dF = Fq{Eo/ E^-'^dE. The flux observed on Earth coming from two images 
formed at a caustic at energy Eo, in the energy interval between 0.9 Eq and Eq, is 

fEo dF A r^o dF 

F^+b(0.9Eo<^<^o)~2 / dE— = ^ 12A / dE— . (2.2) 

Jo.dEo Jl- E/Eq •'^■^^0 dE 



We have neglected a small 0{yl — E/Eq) correction (the term proportional to Ca + 
Cb)- The flux observed from the two images in this energy interval near the caustic is 
12A times larger than the flux of the principal image of the source in the same energy 
range in the absence of magnification. This is also 2.4^4 times the flux that would 
arrive from the principal image at all energies above Eq if there were no magnetic 
lensing. Detection of an UHECR source in a narrow energy range around a caustic 
may thus be more likely than its detection at any higher energy. It may also be more 
likely to detect the CR source at energies around the caustic than at significantly 
lower energies (say half the energy of the caustic or even less) in cases where both 
the principal image as well as the secondary images are not magnified at energies 
below the caustic. 

Large enhancements of the observed flux in a narrow energy range also occur 
in the interesting case when the source position lies along two caustics at nearby 
energies so that the principal image of the source (the only one visible from Earth at 
the highest energies) has also divergent magnification. Such a situation is exemplified 
in figure ^ for a source located at (£ = 90°, 6 = —10°). A pair of images becomes 
visible at energies below Ei^^ = 29.31 EeV, one of which disappears in a caustic 
along with the principal image at an energy E^Q = 22.31 EeV. The secondary image 
that merges with the principal image is that with opposite parity. 

In appendix B we develop the geometrical interpretation and analytic approxi- 
mation to the magnification of images for two nearby folds. The analog result in the 
gravitational lens case was obtained in [0. As shown in appendix B, and verified 
numerically with high accuracy, between the folds the sum of the magnifications of 
the two images with equal parity coincides with the magnification of the image with 
opposite parity, up to a constant term. The analytic fits in figure ^ to the divergent 
magnification of the principal image (/) at E^ and to the magnification at E^^^ of the 
image that survives at low energies (///), lead to Aj = 0.275, Bj = 0.15, Ci = 0.63; 
^iii = -096, Bfjj = —0.11, Cjii = 0.03. The magnification of the image that is 
visible at energies between Eq and E^^ is fit to fiji = + fijjj — 0.7. 
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Figure 2: Example of a source position that is crossed by two folds at nearby energies. 
The solid lines for images I and III are the analytic fits of the numerical results to eq. ( |2.1| ) 
for the amplification near a caustic. The solid line for image II is iii + fini — 0.7. The 
source is located at {i,b) = (90°, -10°). 



3. Angular distribution of events: clustering around caustics 



The arrival directions of the so far detected CRs in the highest energy range is 
compatible (within the limited statistics available) with an isotropic distribution, 
except for some small angle clustering of events (eight doublets and two triplets 
within a total of 92 events with energies above 40 EeV) |Tl|] . The observed relatively 
uniform distribution does not preclude the possibility that UHECRs originate from 
very few nearby sources, if their trajectories underwent significant energy-dependent 
deflections in their journey to the Earth. Trajectories of UHECRs with E/Z below 
approximately 50 EeV are sensibly deflected in the magnetic field model considered 
in this paper. Several quite separated events may actually originate from the same 
source if CRs have a component which is not light (see paper I and references therein). 
It has even been speculated |]12[ that all the events so far detected at energies above 
10^° eV may come from M87 in the Virgo cluster, if the Galaxy has a rather strong 
and extended magnetic wind. 

In paper I we illustrated the angular displacement in the arrival directions of 
CRs from the principal and secondary images of a magnetically lensed source. Here 
we analyse in more detail the angular displacement of image pairs around caustics. 
We argue that caustics can produce a significant clustering of arrival directions of 
events with comparable energies. 

It may well be the case (depending on the sources flux and composition, and on 
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the exact nature of the galactic magnetic field) that a large fraction of the trajectories 
of the UHECRs so far detected have been significantly bent and do not point to 
their sources. In this scenario clustering of events would be infrequent, at relatively 
low energies due to large defiections and at high energies due to the smaller fiux. 
However, when a source is near a caustic, since the significant enhancement of its 
fiux occurs within a narrow energy range, the CRs arrive from relatively nearby 
directions, leading to an angular concentration of events. 

We illustrate the clustering effect in figure |^. We consider the same examples 
of source locations as in the previous section. We assume that the differential flux 
injected by the source scales as and that the detecting system has the same 

efficiency at all energies within the range considered. The energy range was divided 
in 50 bins of equal detection probability. Figure displays the predicted arrival 
directions of the events detected from each source. 

In the case of M87 (top panel) 11 events out of a total of 50 with E/Z larger 
than 10 EeV (half the energy of the caustic) fall in the narrow energy range between 
the energy of the caustic and just 10% below ('+' signs). 3 of those 11 events are 
in the principal image and 8 in the secondary images. Notice that only 5% of the 
events would fall in the same energy range if the source were not magnified. Notice 
also that just 8 events correspond to the principal image at all energies higher than 
the energy of the caustic. The angular clustering effect for M87 is significant but 
not extremely large (in this magnetic field model), partly due to the fact that the 
principal image is also largely magnified around the energy of the caustic, and partly 
because defiections are quite large. 

The middle panel in figure ^ displays the effect for the source located at (£, h) = 
(320°,— 30°). In this case 9 events fall within 10% of the energy of the caustic in 
the secondary images, separated by no more than 5°. Only 4 events are seen for all 
energies higher than the energy of the caustic, scattered over almost 10°. The events 
at energies below the caustic, down to half its energy, are scattered over more than 
20°. 

The bottom panel in figure ^ corresponds to the source located at (£, b) = 
(90°,— 10°), which presents two nearby folds, as depicted in figure 0. In this case 
nearly half of the events (23) fall within the energy range between the two caustics 
(22.3 EeV < E/Z < 29.3 EeV), while 15 events occur for all higher energies, and 
just 12 events appear at lower energies down to E/Z = 7 EeV. The relatively small 
number of events in the lower energy range is due to the large deamplification of the 
image fiux. An unlensed source would have instead 86% of the events in the lower 
energy range considered, 9% in the higher energy range, and just 5% at the inter- 
mediate energies. Notice that while the arrival directions of the events considered 
are spread over more than 15°, nearly half of the events fall within just 2° around 
(/ = 88°, 6= -11°). 

We stress the fact that six out of the eight doublets in the UHECR data listed 
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Figure 3: Predicted arrival directions of 50 UHECR events from M87 (top) and from the 
sources at {e,b) = (320°, -30°) (middle) and {e,b) = (90°, -10°) (bottom), that illustrate 
the clustering effect due to the existence of caustics in the magnetic lens mapping. 



in table 6 in |11] are such that the energies of the events in a pair differ by less than 
10%. The same happens with two events in one of the triplets. This may be an 
indication that at least a fraction of the observed clustering of events may be due to 
magnetic lensing around caustics. 
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Figure ^ (and figure 6 in paper I) illustrate the angular displacement of the 
images as a function of energy. As discussed in appendix A, the observed angular 
separation between images scales clS ■\J~X ^ where x measures the distance from the 
source to the caustic. Consequently, the angular separation between the pair of 
images created at the caustic crossing scales with energy near a caustic as 

M ^ ey^l - E/Eo . (3.1) 

The proportionality constant G is fitted from the numerical output for the angular 
displacement of the images. In the case of M87, O ~ 56°, and in the case of the 
source at {i,b) = (320°, -20°) it is 6 ^ 15.1°. 

4. Effects upon composition 

Another interesting effect of magnetic lensing is that it can strongly alter the mea- 
sured composition of CRs at high energies. This is because at each source location 
the magnification is a function of E/Z. Thus, for a given energy the CR flux of 
components with different Z suffers different (de)magnification as the CRs travel 
through the Galactic magnetic field. 




E (EeV) 

Figure 4: Illustration of the effect of magnetic lensing upon the UHECR observed com- 
position. The source is located at {l,b) = (90°,— 10°), and undergoes the magnification 
depicted in figure |2|. The source injects 50% of protons and 50% of Nitrogen. In the absence 
of lensing the mean value of Z would be {Z) = 4, independently of energy. 

To give an idea of the strength of this effect, let us consider a simple example: 
suppose that from a given source the CR flux arriving to the halo is composed by a 
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fraction / of protons and a fraction 1 — / of heavier nuclei with charge Z^. In the 
absence of magnetic lensing the mean value of Z of the arriving particles would be 
{Z) = /+(! — Due to the effect of the magnetic field, the flux of protons 

on Earth will be (de)magnified by /i(-E') and the flux of heavier nuclei by ^{E/Z). 
Thus, the mean value of Z of the CRs arriving on Earth will be 

^ _ ff^jE) + (1 - f)Z,^^iE/Z,) 

^ ^ ff,{E) + (1 - f)f,{E/Z,) ■ ^ • ^ 

As the magnifications can be large, (Z) can be strongly modified by the magnetic 
field. The effect is specially important for regions of the sky for which the total 
magnification /i (adding the contribution from all the images) has noticeable changes 
as E/Z varies. 

We show as an example the source located at (£, b) = (90°, —10°), which under- 
goes the magnification depicted in figure |^. We took for reference a flux composed 
by a fraction / = 0.5 of protons and a fraction 0.5 of Nitrogen {Zh = 7). In the 
absence of lensing the mean value of Z would be (Z) = 4. We see in figure ^ that 
a clear change from a light composition in the smaller energy region to a heavier 
composition at the larger energies appears due to magnetic lensing. 



5. Time delays 

A charged UHECR that traverses a distance L within a homogeneous magnetic 
field B perpendicular to its trajectory is deflected by an angle of the order of 
rj ~ 5° (10 EeV Z / E){B / ^G){L /kpc) , in the limit of small deflections. Consider two 
initially parallel CRs emitted simultaneously, one with a much lower ratio E/Z than 
the other. If they enter a region permeated by a homogeneous magnetic field and con- 
verge to the same point after a distance L, the CR with lower ratio E/Z arrives later, 
with a relative time delay 6t ^ r]^L/2 ^ 10 yrs(10 EeV Z / Ef{B / ^iGf{L/ kpc)3. 
Higher energy events arrive earlier. 

Time delays induced by intergalactic magnetic fields are an essential ingredient in 
bursting models for the origin of UHECRs . There are very definite observational 
signatures of such a scenario. For instance, individual bursting CR sources would 
have very narrow observed spectra, since only CRs with a fixed time delay would 
be observed at any given time. The energy at the peak in the differential CR flux 
received on Earth should shift with time as 

Here we analyse the implications of the regular component of the galactic mag- 
netic field upon time delays between events from a UHECR source at different ener- 
gies, and between events from different images of a single magnetically lensed source. 

Since UHECRs are extremely relativistic, their time delay compared to straight 
propagation at the speed of light from the source to the observer is simply determined 
by the excess path length. We assume that the extragalactic UHECR sources are 
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sufficiently distant that we can approximate the CR flux incident upon the galactic 
halo as a beam of parallel trajectories, all emitted simultaneously. We then numeri- 
cally determine the difference in path length between the trajectory of a charged CR 
and the path length along the parallel straight trajectory in the beam that reaches 
the Earth. In other words, we measure the delay in the arrival time of a charged 
CR with respect to the arrival time of photons (or charged CRs with much higher 
energy) emitted simultaneously from a distant source. 

Figure § displays the time delay with respect to straight propagation for the 
principal and secondary images of the source locations discussed in the previous 
sections. 

The time delay in normal regions (far from critical points) monotonically in- 
creases with decreasing energy. The spiraled galactic magnetic field is far from ho- 
mogeneous, and thus departures from the E'"^ dependence are certainly expected. 
We have checked that the time delay averaged over a regular grid of arrival direc- 
tions scales, in the BSS-S galactic magnetic field model considered, as < 6t >^ 
1000 yrs(10 EeV Z/EY from very high energies down to E/Z of order 5 EeV. At 
lower E/Z values the increase in the time delay with decreasing energy is faster, 
since a sharp transition from quasirectilinear to drift motion occurs at values of E/Z 
between 3 and 1 EeV. 

What strikes the eye in figure ^ is that the time delay of one of the images in a 
pair can have an energy dependence opposite to that in normal regions. Indeed, the 
time delay of one of the images that are visible at energies below the energy of the 
caustic increases with decreasing energy, while the time delay of the other member 
of the pair decreases. Thus, the relative arrival time of events from a single image of 
a CR source does not necessarily increase with decreasing energy. It is often argued 
that the doublets in which the highest energy event arrived later than the other 
member in the pair can not arise from bursting sources. As we have seen, this is not 
necessarily true near a caustic. 

The relative arrival time delay between equal energy events from different images 
A and S in a pair behaves near the caustic as 



This can be understood as follows, in the limit of small deflections. The trajectory 
of each of the images in the pair is deflected by rj and t] + A6 respectively, with A6 
given by eq. (|3.1|) . The time delays with respect to straight propagation are thus 
proportional to rj'^ and rj"^ + 27] A6 respectively, as long as A6 -C f]. Thus the relative 
time delay between events at the same energy from the two images scales as A6 at 
energies sufficiently close to that of the caustic. The ffi of eq. ( ^.1| ) to the numerical 
output is highly accurate, with T ^ 3200 yrs in the case of M87 and T ^ 1030 yrs 
in the case of the source in the southern sky at [i = 320°, b = —30°). 




(5.1) 
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Figure 5: Time delays in the arrival of UHECRs from the principal and secondary images 
a distant source, measured with respect to the arrival time of a photon (or a much more 
energetic charged CR) . The relative arrival time delay of a pair of images near the caustic 
scales as St oc ^/l — E/Eq. Notice that the arrival time of one of the images decreases for 
decreasing E below the energy caustic, opposite to the behaviour in normal regions. 

6. Conclusions 



UHECRs beyond the ankle in the spectrum are most probably nuclei of extragalactic 
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origin. In their way to the Earth they feel the magnetic field structure permeating 
the Galaxy and hence their trajectories become deflected and their fluxes are lensed. 
The implications of this magnetic lensing are manyfold and they have to be taken 
into account in the analysis of the observations. 

Since the lensing effect depends on the energy it can sizeably affect the observed 
spectrum of the sources. Moreover, for sources located in a large fraction of the 
sky CRs can reach the Earth following different paths, and hence multiple images 
of those sources will appear. A useful way to visualise these effects is to display the 
mapping of the arrival directions in the Earth into the incoming directions outside 
the Galaxy, as was introduced in paper I with what was dubbed a 'sky sheet'. In the 
locations where this surface develops folds, pairs of additional images of the source 
are present. These folds, corresponding to the caustics, move with energy, and as they 
cross a given source location pairs of additional images appear or disappear. Near 



the fold the magniflcation of each image in the pair diverges as /x oc 1/ ^\E — Eq 
Thus the probability to detect events from a given source is noticeably enhanced for 
energies close to £"0, at which the caustic crosses the source location. This also leads 
to an expected concentration of events near the location at which the new pair of 
images appears. This is relevant in the analysis of the small scale clustering present 
in the UHECR distribution and in this respect it is remarkable that the observed 
events in doublets and triplets in most cases are very close in energies. With the 
increased statistics expected with the new detectors, such as Auger [|1^] or High Res 
[H], these features may become testable through a careful analysis of the clustering 



of the events. 

Also the observed CR composition can be affected by magnetic lensing due to 
the dependence of the flux amplification on E/Z. Nuclei with different charges are 
magnified by different amounts for a given energy. This effect is sizeable for sources 
whose magniflcation has a strong energy dependence (in particular when there are 
caustics) and which have a mixed composition. 

Another feature which we have discussed is the time delay due to the galactic 
fleld between the different images of a lensed source. These are typically larger than 
the lifetimes of the CR observatories, and are strongly dependent on the energy. 
If burst sources exist, narrow spectra will then result at any given time, and the 
different images will be observed simultaneously with different E/Z values. 

We have to stress that in addition to the effects related to the magnetic fleld of 
the Galaxy discussed in this paper there may be also similar effects associated to the 
magnetic fleld in the source galaxy or even to the intergalactic flelds if these ones 
are strong. Also the magnetic fleld model adopted here is plausible but the real one 
may differ from it, changing the quantitative details but not the general qualitative 
results. 
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A. Magnification near caustics 

The relation between the source position and the image positions is given by a 
mapping from the source coordinates {£,b)H to the observer's ones {£,b)E- When 
this mapping is multiple valued there will be additional image pairs appearing (see 
paper I). The magnification of a given image will he /i — dQij/dJle, i.e. the ratio of 
differential solid angles associated to the mapping. 

To visualise these facts it is useful to consider the inverse map {£,b)E {£,b)H 
and look at it as the mapping of a surface (the observer's sky) into another surface 
(dubbed the 'sky sheet' in 1), which will be folded when multiple images appear, and 
whose stretching is related to the magnification of the images. The location of the 
folds correspond to the caustics along which image pairs with divergent magnification 
appear and the magnification varies very rapidly with the distance to the fold. Since 
the folds move as the energy is decreased, the knowledge of the magnification as a 
function of the angular distance to the fold for a given energy can be used to obtain 
the magnification, for a given source, as a function of the energy near a caustic, 
which is the quantity of interest to us here. 

The angular dependence of the magnification near a fold has a very simple ge- 
ometrical interpretation in terms of the folded sky sheet. Let us take local angular 
coordinates {x, y) in the source sky such that a; = describes locally the location 
of the fold (and hence the fold is along the y axis while the x axis is orthogonal 
to it) . We can also adopt (non-orthogonal) local coordinates on Earth [X, Y) such 
that X = is mapped into x — Q and similarly y = is mapped into y — Fur- 
thermore, we will assign a third coordinate to the mapping, i.e. {X^Y) — > {x,y,z), 
giving a depth to the sky sheet so that the fold can be visualised. The simplest 
choice for this (arbitrary) third coordinate is to take z = X, as we will do in the 
following. The results will be easier to obtain with this choice, but are independent 
of it. In this way the mapping for a fixed value of F in a neighbourhood of X = 
will look as shown in figure 6. The magnification is given by fx — \KdX/dx\, where 
K=[cosbH\d{e,b)H/d{x,y)\]/[cosbE\d{i,b)E/d{X,Y)\]dy/dY^Ko{Y)+K,{Y)X+ 
0{X'^). The next step is to relate the factor dX/dx — dz/dx with the slope of the 
curve in figure 6, and approximate the fold near the caustic by its Taylor expansion 
X — az'^ + bz^ + 0{z^). Hence we have 

^c^2az + 3bz^ (A.l) 
dz 



14 



and hence, since z = ±^Jx/a{l + 0{x)), with the sign indicating to which side of the 
fold the image belongs, we have 

IJ= ^±B' + 0{y^), (A.2) 

A' and B' being x-independent parameters related to the parameters a, b, Kq and 



Z 




Figure 6: Mapping from the observer's sky {X) into the source sky (x), for fixed Y, with 
the vertical coordinate {z = X) allowing to visualise the fold in the sky-sheet. 

The final step is to relate this with the magnification for one given source as a 
function of energy. Suppose that at energy Eq a source located at (^o, ^o)h would be 
just on top of the fold. If we define the angular 'velocity' of the fold as the energy is 
changed as V{E) = dxf/dE, with Xf being the distance from the source to the fold 
{xf = for E = Eq), one has Xf ~ V{Eq)\Eq — E\. Hence we finally get 

f^{E) = , ±B + 0(^1- E/Eo). (A.3) 

^i^eJeo 

B. Magnification for two nearby folds 

The previous Appendix dealt with the magnification for the two images which appear 
when a fold crosses the source location. The third image (e.g. the one present 
originally) was assumed to be far from the fold so that its magnification was supposed 
to have a non-singular behaviour. Another situation of interest is when two folds are 
nearby (e.g. when the sky-sheet has a narrow fold which moves with energy across 
the source, rather than becoming wider and remaining on top of the source, or when 
the source is near a cusp where two folds merge) so that the three images can have 
large magnifications. 

When a couple of nearby folds crosses a source, we see that a pair of images 
appear at some energy and at a somewhat smaller energy one of the new images 
merges with the original one present at high energies and they disappear (figure |^). 
In this case there is a relation between the magnification of the three images, as we 
now show. A similar relation holds in the gravitational lens case 
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Following similar lines as in the previous Appendix, the folded surface can be 
described as in figure 7. Without loss of generality we can choose the origin x = 
so that the vertical coordinates of the other two images at zi and Z2 in figure ^ 
satisfy Z2 = —z^. In this case a Taylor expansion of the fold will be x = az^ — cz 
(with a, c > and no quadratic term due to our choice of origin for x). The location 
of the folds are at xj = — 2a(c/3a)^/^ and xjj = —xj. The images in region I 

—zj) will have positive parities, while the one 



{z > zj = \Jc/'ia) and III {z < zji = 
in region II {zji < z < zj) will have negative parity. 




Figure 7: Sky-sheet in the source sky corresponding to two nearby folds (for Y fixed). 

The vertical coordinates of the three images for a given x are simply obtained as 

a + 2kTT 



Zk = 2zi cos ■ 



(A; = 0,1, 2) 



(B.l) 



where cos a = x/xjj. Here k = 2 corresponds to the negative parity image (region 
II) while k = and 1 are the images in region I and III respectively. 
The magnification of the images will be /ifc = K\dzk/dx\, where 

dzk 2zi . rcos~^(x/x7/) 

= -^^^=^^^= sin 

da; SxiiJl - {x/xjiY 



+ 



2k7T 



(B.2) 



From this we find that 



dzo 


+ 


dzi 




dZ2 






dx 


dx 




dx 



(B.3) 



i.e. that the sum of the magnifications of the two positive parity images coincides 
with the magnification of the negative parity image. 

Next to leading terms (i.e. corrections of 0{y/x) will add constant terms to the 
magnifications, so that the final result up to corrections of 0{x) is 



f^i + l^iii = fJ-ii + const. 
This theorem is illustrated in figure 0. 



(B.4) 
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